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Abstract 

The non-linear relations between polarization strength and electric field strength for 
ferroelectrics, as well as magnetization strength and magnetic field strength for ferromag- 
netics, can be achieved by introducing the retarded electromagnetic interactions. The 
electric and magnetic hysteretic loops in medium frequency alternating electromagnetic 
fields can be described well. The basic phenomenological formula in current non-linear 
optics can also be obtained in the high frequency fields. The results show that all of these 
non-linear phenomena come from the retarded electromagnetic interaction. 

PACS numbers: 75., 75.50-y, 77.80.-e, 42.65-k 

According to classical electromagnetic theory, the Maxwell equations in medium are 

V-D = pf V-B = (2) 

The experiments show that for general non-ferroelectrics and non-ferromagnetics in weak 
electromagnetic fields, the relations between electric polarization strength P and E, as well as 
magnetization strength Af and H are linear, connected by the so-called constructive equations 

D = eoeE = eoE + P B = pof^H = fio{H + M) (3) 

P = eoXeE M = floXraH (4) 

Here Xe is polarizability and Xm is magnetic susceptibility individually with e = 1 + Xe, 
fi = 1 + Xm- For non-isotropy electric and magnetic mediums, the relations are still linear 
with simple forms 

P = eo^^:E M = po^^ : H (5) 

Here Xe is polarizability tensor and Xm is magnetic susceptibility tensor. 

However, the experiments also show that for ferroelectrics and ferromagnetics, as well 
as in the strong fields, the relations between P and E, as well as M and H are non-linear 
with great complexity. Their relations cant be deduced from theory at present. We now get 
them only by phenomenological hypothesis. In the polarization and magnetization processes 
of ferroelectrics and ferromagnetics, there exist hysteretic phenomena. In the alternating 
electromagnetic fields of low frequency (for example O.lHz), the magnetic hysteretic loops 
and magnetization curve are shown in Figure 1. For ferromagnetics, M is much big than 



H, so the shape of M ~ i7 curve is similar to the shape o{ B ^ H curve. The magnetic 
hysteretic loops are near rectangle. In the figure, the maximum values Mq and Hq are at the 
same point. By fixing the frequency of magnetic field and changing iJo, we can get a series 
of loops. Connecting the vertex of each loop, we get a curve, i.e., magnetization curve. 

When the frequency becomes high, the hysteretic loops become smooth ellipses. When 
the frequency becomes higher, the ellipses become near circular. Figures 2 and 3 are the 
magnetic hysteretic loops of nickel-zinc ferrite in the magnetic fields with frequencies 1500Hz 
and 500Hz individually. 

When the frequencies are high enough, the maximum values Bq and Hq are not at the 
same point. So we cant get the magnetization curves by the same method shown in Figure 

I. By considering continuity of magnetization processes, in medium frequency fields, we 
can estimate the basic sharps of magnetization curves below when the frequency and the 
maximum value of magnetic field are fixed. Because magnetization curves must reach the 
maximum value of magnetic field but cant pass over the value, the magnetization curves can 
take three forms shown in Figure 3, described in broken lines a, b and c. The curves a and b 
are unclosed at first several cycles. It needs some time for them to become stationary loops. 
Curve c can smoothly enter stationary loop at its first cycle before it reaches the maximum 
point of magnetic field. In fact, the situation is the same for magnetization curve in low 
frequency fields. The Loops shown in Figure 1 are only stationary ones. At the beginning 
of magnetization processes, they should also be unclosed. The forms of hysteretic loops and 
polarization curves ferroelectrics are the same as ferromagnetics. 
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Figure 4: 



It is still an unsolved problem how to deduce the function forms of hysteretic loops and 
polarization curves in theory now. What we have at present are the experimental formulas 
for ferromagnetics provided by L.Rayleigh in 1887 

B = t^oifJ'H + bH'') B = fio{fi + bHo)H±,^o^{H^-H^) (6) 

The formulas are only suitable for weak fields. In the formula, fi is initial magnetic suscepti- 
bilities, 6 is a constant. The positive sign is for the up-line of formula of magnetic hysteretic 
loop and the negative sign for down-line. Suppose alternating magnetic field is H = Hq smut, 
let a = hH^/ ^, Eq.(6) can be written as 

B = iJ,oiJ,Ho{smujt + asm^Lot) B = /xo/xi7o[(l + a) smut ± ^cos^cvt] (7) 

Take a = 1, let Hq as the unit of abscissa, 5o/(2mom) ^ the unit of ordinate, we can get the 
magnetic hysteretic loop and magnetization curve shown in Fig.4. It can be seen that the 
loop is not an approximate rectangle, so it dose not consistent with Figure 1 in low frequency 
condition. It is not smooth and continuous at the maximum points, so it dose not consistent 
with Figure 2 and 3 in medium frequency fields. Meanwhile,. Bq and Hq are at the same 
point in figure 3, but the real situation is that they should not for the existence of phase 
difference between them. 

It is proved below that after retarded electromagnetic interaction is introduced, the non- 
linear relations between polarization strength and electric field strength, as well as magneti- 
zation strength and magnetic field strength can be obtained. Based on soft modal theory, the 
electric and magnetic hysteretic loops, as well as polarization and magnetization curves can 
be described well in the alternating electromagnetic fields of low and medium frequencies. 
The phenomenological formula used in non-linear optics can also be obtained by the same 
method in high frequency fields, meaning that the retarded electromagnetic interaction is the 
origin of all these non-linear phenomena. 

The polarization of ferroelectrics is discussed at first. According to the current theory, 
there exist small electrical domain structures in ferroelectrics under critical temperature. 
When the electrical polar moments in domains arrange in parallel, energy is lowest so that 
the system is most stable. Owing to the random distribution of polarization directions of 
electrical domains, when there exist no external fields, whole ferroelectrics arc electrically 
neutral. When external electrical field is added, the polarization directions of domains are 
turned to the direction of electrical field, so that ferroelectrics show strong polarization. In 
alternating electrical field E = Eq sin out, the motion equation of an electrical polar moment 
is 

d^i" ^^dr 9 qEn . 

^ + 2/3-+.oV = ^sm.t (8) 

The second item on the left side of equation represents the damping forces caused by collision 
and radiation and so on with P > 0. The third item represents restoring force and, uq is the 
nature frequency of dipole and q is dipoles charge. The solution of Eq.(8) is 



r = A sin(a;i - ^o) + Ce"^* sin( y'w^ - pH - ■&) (9) 
A= , eo = tg-\4^) (10) 



On the other hand, according to the soft modal theory of ferroelectrics under critical 
temperature, the vibrations of lattices are frozen for the Coulombian forces and restoring 



forces of irons are just offset each other, when there exist no external electric fields. In this 
case, the dipoles would not vibrate with ujq = 0. 

So let cjo = in Eq.(8), the solution of the i-th dipoles motion equation becomes 

n = A sH^t - e,o) + b,- ^e-^^^* A, = — , ° , (ii) 

Here Bi and Cj are constants depending on initial condition. According to the theory of soft 
modal, when t = we have = rjo and Vi = Vio = 0. For simplification, let = 0, by 
relation, tg9io = —2Pi/u, we get 

rio = -^iSm^jO A.i = — Q = - cos 6'io (12) 

So when t < the dipole has a stationery electrical polar moment Pjo = ^i'^io- If ^io is 
positive, 0jo is in the third or fourth quadrant. From tgQiQ = — 2/3j/a;, 9io is in the second or 
the fourth quadrant. So 6iQ should be in the fourth quadrant. Similarly, if rjo is negative, 9io 
should be in the second quadrant. From Eq.(ll) and (12) we get 

Vi = Ai sm{ijjt - Oio) - ^»^Q^^»0 g-2At ^. ^ ^.^ cos((ji - 6*10) - Aiuj cos (13) 

sm UiQ 

The retarded electromagnetic interaction is introduced below. Let t, fi, t/j, a, and a, rep- 
resent retarded time, coordinate, velocity, acceleration and acceleration of acceleration, t, fi, 
Vi, Si and Si represent non-retarded time, coordinate, velocity, acceleration and acceleration 
of acceleration. A particle with charge qi , velocity Vi and acceleration Si at space point fi (t) 
and time t would cause retarded potentials at point fi{t) and time t as follows 

47r£o(l-^)n, 47r6oc2(l - ^)r., 

In the formula fij{t,t) = fi{t) — fi{t), t = t — rij{t)/c, Vij =| fij \, fij = fij/rij. When 
Vi << c, the retarded distance rij{t) at retarded time t can be replaced approximately by 
non-retarded distance rij{t), i.e., we can let t = t — rij{t)/c t — rij{t)/c and get 

rij{t, t) =1 fj{t) - fi{t) 1=1 fj{t) -n[t- rij{t)/c] HI fj{t) -fi[t- rij{t)/c] \= njit) (15) 

It is noted that ^^(t) 7^ Tij{t), for r,y(t) is the approximate retarded distance, but rjj(t) is 
not the retarded distance. In this case, wc can develop retarded quantities into series in light 
of small quantity r^j/c. By relation Vi = dfi/dt, Si = cPfi/dt"^, Si = d^fi/dt^ ■ ••, we get 



fij{t, t) ~ fi{t) - fj{t - nj/c) = fi{t) - fj{t) + ^n,{t) - ^r^it) + ^rf^it) + 



nAt) + '-fnAt) - '-^r^t) + '-§-rUt) + • ■ ■ (16) 



Vi{t) vi{t) - ^njit) + ^4(t) + ■ • • Si{t) it) - ^Tijit) + ... (17) 
Thus, according to Eq.(16), the retarded electrical dipole moment becomes approximately 



When t is big enough, we have e ^'^^^ 0, the second item on the left side of Eq.(13) can be 
omitted, so we have rij = Ai sin(wt — ^o) ^ well as 

Vi = AiU cos{ut — Oio) CLi = —AiU^ sin(a;i — Oio) = —Aico^ cos{ut — ^jo) (19) 

If there is only one electrical dipole, the directions of fij, Vi, Si and are on the same straight 
line. In this case, Eq.(18) can be written as 

Pi = — 2 [sin(u;t - 6*10) + aj-Eo cos(cjt - 9io) | sin(wi - 9io) | 

sin(a;t - ^^o) I sin(u;t - ^^o) P ^ cos(wt - 9io) \ sm{u;t - 9io) f + • ••] (20) 

Here ctj = Qi cos 9io/{cmiLo). In this way, the non- linear relation between polarization strength 
of ferroelectrics and electrical field strength is introduced automatically. The order of mag- 
nitude of non- linear items or the value of parameter = qi/{cmiuj) is estimated below. 
For ferroelectrics composed of BaTiO^, mi = 8.02 x lO^^^Kg is the mass of Tj iron with 
Qi = 4e. Let 9io = for simplification, we can get Ui = 2.66 X 10-V'^- Taking a; = 27r x 10^, 
Eo = 10^ V/M, we have aiEg = 0.42. The non-linear effect is quite strong. Similar, for O 
iron, mi = 2.67 x 1.67 x W^^^Kg, qi = — 2e, we have UiEg = —0.64. If frequency becomes 
lower, the non-linear eS"ect would become bigger, for example, taking lo = 207r, Eq = W'^V/A, 
we have aiEo = 4.2 for Ti iron and ctj-Eo = 64 for O iron. When frequency becomes high, 
the non-linear items become small, when frequency become high enough, the non-linear ef- 
fect would become small enough so that it cant be observed. The frequency is just so-called 
cut-off frequency. On the other hand, if the natural frequency of dipole loq 7^ 0, we have 
coq = 10^^ ~ 10^^ >> CO. When electrical field is not very strong, for example Eq = W^V/M, 
we have aiEo = 10~^ ~ 10~^, the non- linear effect is too weak to be observed. So the 
hypotheses of lattice vibrating frozen in soft modal theory is important for us to get observ- 
able non-linear relation. It is obvious that all results above coincide with experiments. The 
parameters of each item in Eq.(20) are ib(aj-Eo)"'/n!. Because aiEo is finite, when n ^ 00 
we have {aiEoy-fnl 0. Write the n item of the series as Sn, we have Sn+i/Sn ~ aiEo/n. 
When n ^ 00 we have Sn+i/Sn 0, So the series is convergent. Though when frequency 
is low and electric field is not very high and n is not very big, wc can have (aiEo)^ /n\ > 1, 
even {aiEo)"^ /n\ » 1. In this case many items are needed to engage enough accuracy. 

Eq.(20) is based on simple model, but practical ferroelectrics are very complex. For 
practical ferroelectrics, many factors should be considered, for example, the practical lattice 
structure, pseudo spin waves and statistical nature of system, interactions between dipoles, 
as well as interactions between irons and electrons and so on. We can put all of those factors 
in a nutshell as a factor ki, and write the polarization relation of practical ferroelectrics as 

Pi = — :r-^[kii sm{ut - 9io) + k2iaiEo cos(a;t - 9io) \ sm{ujt - 9io) \ 

rriiU^ 

^kuo^ sin(a;i - ^^o) I sin(u;t - 0,o) |' cos(a;t - 0,o) I sm{ut - 9io) \^ + ■ ••] (21) 

Suppose there are N effective dipoles in ferroelectrics, the total polarization strength of 
ferroelectrics in external alternating electrical field is 

P = 9 ' [kii sm.{ut - 6*10) + k2iaiEQ cos(wt - 9iQ) \ sm{iot - 9iQ) \ 

rriiOJ'^ 

^hi^ sin3(a;t - ^,0) - cos(a;t - ^,0) I sm{u;t - 9io) f + • ••] (22) 



The first item can be written as 



— 5 sin(u;r — t^jo) = £o\Xi svnut — yi cos u)t) [26) 



■ 1 rriiOJ 
1=1 * 



2 — yi = 2 (2^) 

In the formula, xi and yi can be regarded independent each other, so we regard them as the 



ordinates of the points in a plane. Let Xei = yXi+yf, yi/xi = tgOi, we can always write 
xi and yi as xi = Xei cos^i and yi = Xei sin^i. So Eq.(23) can be re- written as 

qfEokii cos Oio . , . q . • / + a \ /r,K^ 
2 sm(a;t - Oio) = £oXei sm{Lot - Oi) (25) 

Similarly, the second item in Eq.(22) can be written as 

^q^k2iaiEl cos Oio a \ \ ■ i ^ a 

> Ti cos(a;i — Pio) sm(a;i — On 

^ — ^ Tin ■/, 

1=1 



= £o {x2 sin Lot COS Lot — 2/2 COS ut sin oot + sin^ ojt — Z2 s\v? ut) (26) 
A qfk2iaiE^ cos^ \ cos 6*10 | -A qfk2iaiE^ sin 6*10 cos 6*10 | sin 6*10 | 

^2 = 2^ 2 y'i = l^ 2 (2 V 

g,^A;2iai^o sin^iocos^io I cos6'io I g-^fejajE'o cos^ 6*^0 | sin 6*^0 | 
2 ^2 = V 2 

Because of sin^jo = cos(7r/2 — ^jo)) it can be considered to be the same to take sum over 
sin^io I cos^io I and over cos^jo I sinSjo I- So we have -22 = -^2- In this way, there are three 
independent parameters on the right side of Eq.(26). We can take them as X2-, and 
write Eq.(26) as 

J2 cos(cc;i - Oio) I sin(a;i - ^^o) |= eoXe2^o cos(u;t - ^2) | sin(u;t - ^3) I (29) 

i=i 



Similarly, the third item has four independent parameters xs, ^4, ^5, 9q, and the third item 
has five independent parameters Xit ^Ti thetag, 9g, 610, and so on. At last, Eq.(22) can be 
written as 

P = eo-Eo[xei sin(a;t - 6*1) + Xe2-Eo cos(a;t - 62) \ s'm{ujt - 63) | 

+Xe3-E'o sin((ji — 64) | sin(a;t — ^5) sin(u;t — 6q) 

+Xe4-£'o cos(a;t — O7) \ sin(a;t — thetas) sin(a;i — ^9) sin(a;t — ^10) | + • ■•] (30) 

The parameters Xen and 0„ are different for different ferroelectrics and can be determined by 
experiments. The formula can also be written as 

P = £oXei-£'o[sin(a;t — ^i) + bi cos{ujt — 62) \ sin(a;t — ^3) | 

+62 sin(ajt — 6,i) \ sm{ujt — 0^) sin(u;t — 0q) \ 
+63 cos(a;t — 67) I sin(wt — 6*8) sin(t(;i — ^9) sin(u;t — ^10) | + • ••] (31) 
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Figure 5: 

Here 6„ = Xen-iEQ~^ /xei with 6„ ^ when n oo. Then, wc consider the simplest 
situation to take two items in (31). Let bi = 0.1, bn = 0(n > 1) and Oi = 02 = Oz = When 
Lot = 9 we get P = 0. In this case 9 is just the polarization hysteretic phase. Take 9 = 30f 
and 60^^. Tate Eq as the unit of abscissa, let Pq represent the maximum value of P and take 
Po/i^oXei) as the unit of ordinate, wc get the hysteretic loops shown in Fig. 5 and 6. They 
coincide with practical results in the medium frequency fields. 

For low frequencies fields, we stochastically take bi = 0.400, 62 = 0.256, 63 = 0.461, 
64 = 0.164, 65 = 0.665, 66 = 0.184, 67 = -0.419, bg = 0.0420, 69 = 0.034, 610 = 0.027, 
b^y^Q = 0, 6*1 = 0.300, 6*2 = 2.100, 6*3 = 1.500, 6*4 = -0.300, 6*5 = 1.300, 6*6 = 6.000, 
6»7 = 6.000, 6*8 = 1.200, 9n>s = 1-050 + {2n)mod6, get Fig.7. Similarly, take bi = 0.464, 
62 = 0.147, 63 = 0.067, 64 = 0.164, 65 = 0.475, be = 0.021, 67 = -0.042, bg = 0.042, 
69 = 0.034, 610 = 0.027, 6„>io = 0, ^1 = 0.270, 92 = 3.180, ^3 = 1-780, ^4 = 0.080, 
6I5 = 0.160, 6*6 = 3.030, 97 = 5.580, 6*8 = 5.560, 9n>8 = 1 + 0.200n, get Fig.8. Two of them 
are good rectangles, showing that Eq.(31) can also be used to describe the hysteretic loops 
in low frequency fields. 





Figure 7: 




Figure 8: 



As for polarization curves, by using Eq.(12), Eq.(22) can be written as 



qfkuEo cos 9io . . , . ^ cos,^ 



^ mnjj^ sm 

g|te£|cos^ - ^,o) - cos0,oe-^^^*] | sin(a;t - 0,o) - ^e'^f'^' \ • • ■ (32) 

Similar to Eq.(30), by considering relation 2/3jt =| tgQi^ \ u)t, the formula above can be written 
as 

P = eoXeiEo{sm{ut - 9^) + g^e-'''\t90i\->t 

+bi[cOs{ut - 02) + ^26-^=21*^^^1'^*] I sij^(^i _ + g^e-kz\tge2\o:t || . (33) 

Though each electrical domain is polarized, whole ferroelectrics are medium generally at 
beginning. The initial condition is P = when t = 0. Only taking first two items, we have 
gi = sin 9i, gz = sin ^3 and get at last 

P = eoXeiEo{sm{ut - 9i) + sin0ie-'=il*s^il'^* 

+6i[cos(wt - 62) + (72e-'=2l*»^2l'^*] I _ Q^-^ + sin03e-'=3|t9e2|a;t II (34) 

The parameters ki > 0, the value should be decided by experiments. Taking Oi = O2 = 9^ = 
9 = 3GP, b = 0.1, ki = k2 = ks = 1, 32 = 1, we get polarization curve shown in Fig.5. Taking 
9 = 60^, b = 0.1, ki = k2 = ks = 0.1, g2 = — cos 60^^ = —0.5, we get polarization curve shown 
in Fig.6. 

The results are the same for the magnetic hysteretic loops and magnetization curves of 
ferromagnetics. As we known that magnetic phenomena can be equivalently described by 
the concept of magnetic charges. According to quantum mechanics, there exist the exchange 
forces between charges. When magnetic moments in magnetic domains are arranged in 
parallel, the systems energy is lowest so that the system is most stable. Magnetic moments are 
caused when electrons move around atomic nucleuses. If there is no external field, electrons 
motions are stationary and the magnetic moments of atoms are also stationary without 
varying with time. The situation is similar to the soft modal theory of ferroelectrics. So we 
can thinks that the vibration of magnetic dipole is frozen with nature frequency uiq = when 
external fields do not exist. In this way, Eq.(ll) can also be used to described the motion of 
magnetic dipole in alternating magnetic field. Similar to Eq.(20), in the alternating magnetic 
field H = Hq smut, the magnetic moment of a single magnetic dipole is 

Mi = — ^3 [sin(a;i - 9io) + UiHo cosiut - 9^) | sin(wi - 6*10) I 

+^p- sm{ujt - 9io) I sin(u;t - ^^0) P ^ cos{ut - 9io) \ sm{ujt - 9io) \^ + ■ ••] (35) 

The order of magnitude of parameter ai = qi/{cmiuj) is estimated below. It is known that 
only out-shell electrons have contribution to magnetic moments. Suppose there is only one 
out-shell electron in an atom, each atom has a magnetic moment called Bohr magneton 
Mi = ^Qeh/{2ms) ~ 1.2 x lO^^^Wb ■ M. Suppose the distance between atoms is about 
~ 10~^°iV ■ M/A, according to the definition of we have 

the value of equivalent magnetic charge qi ~ 1.2 x lO^^^M^. The value is similar to electrons 
charge. If atomic nucleus or iron is considered at static and electron moves, rrii should be 
electrons mass, we have aj = 4.4 x lO^/w, much big than that of ferroelectrics. Taking 
a; = 27r X 10^, Hq = O.IA/M, we have aji^o = 0.7 with strong non-linearity. If rrii is regarded 



as the masses of atomic nucleus or iron, the order of magnitude of is similar to that of 
ferroelectrics. In practical situations, effective masses will be taken. Similar to Eq.(30), we 
have 

M = i^oHolni sm{ujt — 6i) + IJ.2H0 cos{ujt — 62) \ sm{ujt — 63) \ 

+/X3iJo sin(a;t — 64) \ sin(a;t — ^5) sm{ujt — 6e) \ 

+/X4i?o cos{ujt — 6y) I sin(a;i — 6s) sin(a;t — ^g) sin(a;t — ^lo) | + • •■] (36) 

The magnetic hysteretic loops can also be described well by the formula above. We can also 
get similar formula to describe magnetization curves of ferromagnetics, but discuss no any 
more here. 

The same method can be used to describer the non-linear phenomena in optics. In 
the current non-linear optics, interaction between light and materials is described by half- 
classical method, i.e., to use classical theory of electromagnetic field to describe lights motion 
and quantum theory to describe interaction between photons and material particles. The 
non-linear phenomena can be dealt with well by the non-linear polarization formula 

P = eoiXeiE + Xe2 ■ EE + : EEE ■ ■•) (37) 

Because optical crystal is not ferroelectrics in general, we regard atoms and molecular 
as dipoles with nature frequency loq 7^ 0. The motion equation and solution of dipolcs are 
Eq.(8) and (9). When t 00, the polarization can also be described by Eq.(30). Taking 
61 = 03 = 65 = 06 = 07 = 08 = 09 = 010 = 0-, theta2 = 07 = + 37r/2, in spite of the phase 
difference 6 thought it exist actually, we get 

P = eo{XelE + Xe2-EE + Xe3E^E--) (38) 

The result is similar to Eq.(37), from it most of non-linear phenomena just as multiple 
frequencies, sum frequencies and different frequencies and so on will be described. If non- 
isotropy of crystal is considered, we should write (38) as (37). Tt should be noted that 
Eq.(37) is only a special situation of Eq.(30). 

Because nuclear mass is much big than electrons mass, the displacement of electron caused 
by external field is much big than nuclear one, so we can think that nucleus is at rest in 
center and electron vibrates around center. Suppose the wavelength of incident light is 
rrii is electron mass, we have = qi/{cmiUj) ~ 1.2 x 10 ^^M/V The 
strength of electric field of laser used in non-linear optics is about Eq ~ IQ^^V/M in general, 
so aiEQ ~ 10~^means that the radio between the first item and second item in Eq.(38) is about 
10~^. The value coincides with experiments with aj = 10~^^ ~ 10~^^, aiEo = 10~^ - 10"^^^^. 
If the light used is not laser with £^0 < W'^V/M in general, we have UiEq < 10"^, the 
non-linear effects are too weak to be observed. 

In order to explain the origin of non-linear phenomena in optics and get Eq. (37) in theory, 
the non-linear oscillator model is used at present. The motion equation of oscillator is written 
as W 

_ + 2/,-+u,„V + 0.' = - (39) 

In the formula, mDr^ is a non-linear force to cause non-linear effects. However, it can be seen 
that the origin of this item is still unclear. Where it comes from? Why does it not taking the 
form such as mDr^^^ or mDr^ and so on? The current theory cant provide clear explanation. 
Meanwhile, the value of parameter D cant be decided by theory now. Quantum mechan- 
ics provides a method of perturbation approximation to calculate polar izability, but cant 
yet explain the origin of non-linear phenomena. It is obvious that by introducing retarded 
electromagnetic interaction, we can explain the origin of non-linear effects in ferroelectrics. 



ferromagnetics and optics well in a united form from the macro-angle of classical electromag- 
netic theory. The result also prompts us to consider retarded interaction in quantum theory 
to get more accurate description of ferroelectrics, ferromagnetics, non-linear optics and other 
fields. 

On the other hand, it can be seen that Eq.(14) can not keep unchanged under time reversal 
with Vi —Vi, meaning that retarded electromagnetic interaction is asymmetry under time 
reversal. It is the same in Eq.(16) under time reversal with Vi — > —Vi, di — >■ Oj, — > — Oj, 
meaning that r also violates symmetry of time reversal. So in the non-liner processes of 
ferroelectrics, ferromagnetics and non-linear optics, symmetry of time reversal is violated 
obviously. In fact, the polarization and magnetization processes that electromagnetic fields 
change along positive directions are just the time reversal processes in which the fields change 
along negative directions. The hysteresis loops shapes are similar to that in the cycle processes 
of heat engines, both are irreversible to produce dissipative heat. For the polarization and 
magnetization processes of ferroelectrics and ferromagnetics, only electromagnetic interaction 
is involved. So the irreversibility is caused only by electromagnetic interaction itself, or 
speaking strictly, by retarded electromagnetic interaction. It is impossible to explain this 
kind of irreversibility by the theories of coarse graining and mixing current and so on as 
advocated in current statistical mechanics. 
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